Quantum Dot in the Kondo Regime coupled to p-wave superconductors 
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This paper studies the physics of junctions containing superconducting (S) and normal (TV) leads 
weakly coupled to an Anderson impurity in the Kondo regime (K). Special attention is devoted to 
the case where one of the leads is a p — wave superconductor where mid-gap surface states play an 
important role in the tunneling processes and help the formation of Kondo resonance. The novel 
physics in these systems beyond that encountered in quantum dots coupled only to to normal leads 
is that electron transport at finite bias eV in SKN and SKS junctions is governed by Andreev 
reflections. These enable the occurrence of dissipative current even when the bias eV is smaller than 
the superconducting gap A. Using the slave boson mean field approximation the current, shot-noise 
power and Fano factor are calculated as functions of the applied bias voltage in the sub-gap region 
eV < A and found to be strongly dependent on the ratio tx between the Kondo temperature Tk 
and the superconducting gap A. In particular, for large values of tx the attenuation of current due 
to the existence of the superconducting gap is compensated by the Kondo effect. This scenario is 
manifested also in the behavior of the Josephson current as function of temperature. 
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INTRODUCTION 

The Kondo effect is one of the simplest manifestations 
of many-body physics exhibiting strong correlations |jj . 
Originally, it has been introduced in order to explain 
the occurrence of a shallow minimum of the resistivity 
p(T) of bulk metals doped with a small concentration of 
magnetic impurities. It has been confirmed that at low 
temperatures (below the Kondo temperature Tk) scat- 
tering from magnetic impurities is strong, leading to an 
enhanced resistivity. It is now well established that at 
this strong coupling regime the theory of perturbation 
fails due to logarithmic divergences of higher (than sec- 
ond) order diagrams. 

The Kondo physics also plays an important role in elec- 
tron transport through quantum dots, where instead of a 
magnetic impurity one now encounters localized electrons 
H ||, |J . Inspection of the dot-electron Green function re- 
veals an enhancement of the local density of states at the 
Fermi energy (the Abrikosov-Suhl resonance). Thus, un- 
like the case of Kondo effect in bulk metals, its hallmark 
experimental manifestation in quantum dots is an en- 
hancement of the conductance. An important milestone 
in this field has been recorded by recent observations of 
the Kondo effect in transport through quantum dots|||, ||] 
and in a finite length carbon nanotubes [Q . These experi- 
ments have paved the way for a new class of experimental 
investigations into the physics of strongly correlated elec- 
trons in general and the Kondo effect in particular. To- 
gether with more recent experiments || demonstrate the 
feasibility of exploiting tunable physical parameters of a 
quantum dot system in order to yield important informa- 
tion on Kondo systems and other many-body related phe- 
nomena. Such studies, for example, encompass the en- 
tire crossover region between the Kondo limit, the mixed 
valence regime and the non-Kondo (weak-coupling) do- 
main. 

While the Kondo physics in a quantum dot attached 



on both its sides to normal (N) metallic leads received 
much recent attention, it has been realized that novel 
physical effects emerge if (one or both) electrodes at- 
tached through a tunneling barrier to a quantum dot in 
the Kondo regime (K) is a superconductor (S). The 
central electron-transport mechanism in such SKS and 
SKN junctions is that of Andreev reflections when two 
particles tunnel together coherently to form a Cooper 
pair in the superconductor. In fact, the physics of SKS 
and SKN junctions is determined by the interplay of 
these Andreev reflections with the Kondo-formed reso- 
nance in the spectral density of states of the dot elec- 
tron |||. This interplay between the Kondo effect and 
superconductivity in the sub-gap region is rather effec- 
tive here. We then expect the the I — V characteristics 
in the sub-gap region to be most important. Andreev 
reflections play a dominant role both in and out of equi- 
librium. In the former case they are responsible for oc- 
currence of direct Josephson current while in the latter 
case they are the cause of dissipative currents at sub-gap 
voltages as well as for suppressing the zero bias anomaly 
in NS junctions. The effect becomes even more signifi- 
cant when the effective Kondo temperature Tk exceeds 
the superconducting energy gap A. A natural candi- 
date for such SKN junction is that in which the role 
of the quantum dot or the Kondo impurity is played by 
a Carbon nano-tube (CNT), where relatively high val- 
ues of Tk can be achieved. Indeed, Kondo effect with 
rather high Kondo temperature (Tk > IK) has been re- 
ported in N — CNT — N junctions with a few conducting 
channels (7). Moreover, superconducting junctions with 
a weak link formed by CNT have already been fabricated 

MM- 

Let us illuminate the somewhat subtle distinction be- 
tween s — wave and p— wave SKN junctions. In the case 
of an SKN junction when one electrode is an s — wave 
superconductor and the Kondo impurity is weakly cou- 
pled to the S and N electrodes, is determined by the 
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competition between two important effects. The first 
one is the formation of a Kondo singlet which screens 
the bare impurity spin and drives the system toward the 
unitary limit at very low temperatures. The second one 
is the existence of the superconducting gap which im- 
plies a vanishingly small density of low energy electron 
states [[l2| [l| |l4| . These are precisely the electron states 
which are needed in order to screen the Kondo impurity. 
Note that if the Kondo temperature is larger than the 
superconducting gap (Tk > A) then electron states out- 
side the gap can participate in the screening interaction. 
However, this condition cannot be easily realized. Con- 
sider, on the other hand, an SKN junction in which the 
S electrode consists of an unconventional superconductor 
which is oriented relative to the interface in such a way 
that the pair potential reverses its sign on the Fermi sur- 
face. In this case, mid-gap zero-energy states (ZES) are 
formed which are localized near the surface of the uncon- 
ventional superconductor. These states can now partic- 
ipate in screening the impurity spin through the Kondo 
effect and emergence of sub-gap current is expected. In 
the experimentally feasible setup of S — CNT — N where 
just a few tunneling channels are present, charge is car- 
ried mainly by quasiparticles moving perpendicular to 
the interface. This restricts the possible values of the 
angle 9 between the superconducting surface and the di- 
rection of the injected quasiparticles Formation of 
ZES is possible when A(0) = -A(tt - 0). If the impu- 
rity is almost point-like, the relevant injection angle is of 
course 9 — 0. Unconventional superconductors for which 
the symmetry of the pair potential is that of a triplet 
(p — wave superconductors) satisfy the above condition. 
These have recently been discovered by Maeno, et al 
in Sr2RuO<±. For impurities of finite extent, one may 
also consider formation of ZES in d — wave supercon- 
ductors. Thus, the physics of SKN junctions with an S 
electrode whose order parameter has a non-trivial sym- 
metry is affected by the formation of ZES in the Kondo 
regime. 

Beyond investigating the conductance dependence on 
the applied bias we propose and explore theoretically 
other novel experimental tools suitable for probing the 
Kondo regime. These include shot-noise measurements 
and the Josephson (direct) current. Like the conduc- 
tance, we will calculate and analyze shot-noise power 
spectrum in SKN junctions at very low temperature 
where one lead is either an s — wave superconductor 
(mid-gap surface states are absent) or a p — wave super- 
conductor (mid-gap surface states are present). We will 
also study the temperature dependence of the Josephson 
current in these junctions. 

In order to study this strongly interacting non- 
equilibrium problem in the framework of workable ap- 
proximation, we restrict ourselves to voltages that do not 
exceed the superconducting gap so that the Kondo tem- 
perature can satisfy the inequality Tk > max(A,eV). 
Under this condition we can apply the well known ap- 
proach familiar in the physics of the Anderson impu- 



rity Hamiltonian |17j], namely the mean field slave bo- 
son approximation (MFSBA) (see below). It consists of 
a workable scheme for calculating the conductance and 
the zero frequency shot-noise power. (Here we also ap- 
ply the MFSBA to calculate Josephson current in SKS 
junctions.) 

The combined effect of multiple Andreev reflections 
(MAR) and electron-electron interactions on the shot- 
noise in SKS junctions was studied in Ref. g |2C]. 
It was shown that in the strong coupling Kondo regime 
the effective transmission is enhanced, and a ballistic-like 
channel opens up. Hence, an interplay between MAR and 
the Kondo resonance yields an excess current at zero bias 
and the I — V curve behaves similarly to that of nonin- 
ter acting ballistic junctions. At very high values of Tk 
and in the low voltage limit the current approaches the 
nonintcracting value Iar = -ieA/h. Analogously, the 
shot-noise power is shown to display a pronounced maxi- 
mum at V = and a decay as 1/V at small bias, familiar 
in the standard theory of noninteracting SNS junctions 

In section II the Hamiltonian of SKN and SKS junc- 
tions is defined, and the bare action is introduced. After 
integrating out the fermion fields pertaining to the super- 
conducting leads an effective action is obtained, depend- 
ing solely on the dot variables. Calculations and presen- 
tations of conductance, shot-noise power and Josephson 
current are respectively detailed in sections III, IV and 
V. 



MODEL HAMILTONIAN AND EFFECTIVE 
ACTION 

The dynamics of systems like SKN or SKS junctions 
is governed by the Hamiltonian 

H = H L + H R + H d + H t + H c , (1) 

in which Hj, (j = L,R) are the Hamiltonians of the 
electrodes which depend on the electron field operators 
ipaa{^, t) where r = [x, y) and a = ± is the spin index. 

i = / dr[*t <7 (r)e(V)* i -7*} t (r)*} i *ii(r)*i T (r)].(2) 

Here 7 is the BCS coupling constant and £(V) = 
— V 2 /2m — n with /x being the chemical potential at tem- 
perature T. The Planck constant is set H = 1 and when- 
ever appropriate, the spin, space and time dependence of 
all the field operators will not be explicitly displayed. As 
in Refs. |l8|, ^ the dot is represented by a single level 
Anderson impurity with energy eo < and Hubbard re- 
pulsion parameter U. In the Kondo regime of interest 
here we set C/ — ^ 00 and assume | eo | to exceed any other 
energy scale except U. In this case it is convenient to 
express the dot and the tunneling Hamiltonians Hd and 
H t via slave boson (operators 6, w) and slave fermion 
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(operators c, c 1 ) auxiliary fields p3[. Explicitly, 

a 

H t = 5^054^(0, t) + /i.e., (3) 

jo- 

where Tj is the tunneling amplitude. Finally, the Hamil- 
tonian of the system must also include a term which pre- 
vents double occupancy in the limit U — > oo. This term 
reads, 



(4) 



where A is a Lagrange multiplier. 

Following Ref. jl8| let us consider the dynamical "par- 
tition function" 



V[F] exp(iS), 



(5) 



where the path integral is carried out over all fields [F] 
and the action S is obtained by integrating the La- 
grangian pertaining to the Hamiltonian ([!]) along the 
Keldysh contour. In performing the functional integra- 
tions the boson field operators are treated as c-numbers. 
As a result one arrives at an effective action expressed in 
terms of the Green functions of the leads. 

S c ff = -iTrlnG- 1 - / dt[Xa z (bb - 1)]. (6) 



Here A = (Ai,A2), 6 = (61,62) and a z are diagonal ma- 
trices acting in Keldysh space. The inverse propagator 
G _1 depends on the Green functions of the electrodes. 22 
Performing the standard basis rotation in Keldysh space 
one finds, 



G-\e,e')=8{e-e'){e 



T6 2 



T z g+(e)T z ), (7) 



where e = eo + A is the renormalized level position (in the 
Kondo limit one has e ~ 0) and T = {T L +T R )/2 cx T c .n 
is the usual transparency parameter. The 2x2 matrix 
representation (in Keldysh space) for g is composed of 
diagonal elements g R ' A {e) and an upper off-diagonal el- 
ement g K (e) = (g R — g )th(e/2T), Here and below we 
define 



g± = ilHl ± iBgR., 



(8) 



with asymmetry parameters jj = Tj/T. The matrix gn 
has the standard structure with retarded and advanced 
superconductor Green functions which in the s — wave 
case reads, 



(e±iO) + |A|r B 
V(e±*0) 2 -|A| 



(9) 



while for the p — wave case with incidence angle a it 
reads, 



with 



(10) 



R/A 



<h 



-R/A 
92 = 



i i/(e ± 0) 2 — I A| 2 cos a — e sin a 
e cos a + i y/(e ± 0) 2 — | A| 2 sin a 
ilAI 



ecosa + iyj (e ± 0) 2 — |A| 2 sin< 



(11) 



Note that the Pauli matrices T XjVtZ act in Nambu space. 
For a one channel system one can take a=0 and simplify 
the expression for the Green functions, 



.KM 1 
9i = - 



e±0 



-R/A _ »|A| 

92 ~e±0' 



(12) 



In the case of NKS junction, the left lead is represented 
by the Keldysh Green's function of a normal metal g^: 

9l ~ ^ ! l ne Keldysh Green function g^ nas onr y di- 
agonal matrix elements equal to 2ith[(e ± eV)/2T], Per- 
forming the variation of the effective action with respect 
to the fields 6 and A a couple of self-consistency equa- 
tions are obtained that determine these fields. In or- 
der to explicitly write down these self-consistency equa- 
tions let us introduce the bare Kondo temperature T K = 
Dexp[— 7r||eo|/(2r)] and define a parameter X by T6 2 = 
T K X, where D is the energy bandwidth. Then the 
MFSBA equations take the form 



X = — 1 -^—TtG k t z , 
iT 



(13) 



A = -Tr[G K r z (g R + g A ) + (G R + G^t^t^A) 

where the trace includes energy integration as well. 
Eq. (|ll]) effectively determines the Kondo tempera- 
ture (through the parameter X), and reflects the con- 
straint which prevents double occupancy in the limit 
U — ► 00. The second self-consistency equation (|l4| ) de- 
fines the renormalized energy level position e. Let us 
briefly discuss the validity range of the present analysis. 
The MFSBA is known to encode the Kondo Fermi-liquid 
behavior at low temperatures. An important parameter 
here is the ratio between the Kondo temperature and the 
superconducting gap t K = T K /A ^ ||, [[§. For t K > 1 
a Fermi liquid behavior is expected. Accordingly, in this 
regime Eq. ([l3]) should have a nonzero solution 1^0 
which corresponds to nonzero Tk- On the other hand, in 
the limit of large A the only possible solution is the triv- 
ial one 6 = (and, hence, Tk = 0). Quantitatively, the 
MFSBA is reliable only for sufficiently large values of tK- 
We believe, however, that it can provide useful qualita- 
tive information also for moderate values of tK describing 
a crossover between the Kondo regime and the Coulomb 
blockade dominant domain 123] . It is worth noting here 
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that, strictly speaking, the applied bias voltage V attenu- 
ates the Kondo resonance and lowers Tk ■ Hence, for the 
reliability of the MFSBA in non-equilibrium situations, 
both A and eV should not exceed the Kondo tempera- 
ture. Attention below is mainly focused on the sub-gap 
voltage regime eV < A in which case tic appears to be 
the only relevant parameter. 




CONDUCTANCE 1 1 r 

0.0 0.5 1.0 

eV/A 

The expression for the tunneling current obtained in 
Rcf. [fijl has a simple representation in terms of the dot 
Green function, 



I = i^Tr[(G R T z -T z G A )g«-G K g], (15) 

Fig. 1 The conductance G (in units of e 2 /h) versus the bias 
V (in units of A/e) for an s — wave SKN (dash curves) and 
where for SKN junctions we denote Sp - waveKN (solid lines) junction at sub-gap voltages with 

T/Tk = 200. The parameter t K =2,3,5,100 (from down to 
top). The upper line corresponding to ti< =100 coincides for 
g = —"/n{g R T z — T z g A ) — 2i"/LT z , s and p wave superconducting leads. 

~9 R ' A = -Mt z 9 R > A + 9 R > A t z ). (16) 



The last quantity will be used later on as it enters the 
expression for the noise power spectrum. Being com- 
bined with eqs. (H) and @ the result (§J) can be 
conveniently used for computing the transport current 
and the differential conductance of an SKN junction in 
the Kondo regime for different values of tK ■ 

For sufficiently large tK we anticipate a strong Kondo 
resonance and the G — V curves for both s — wave and 
p — wave superconductors is expected to resemble that 
of purely ballistic junctions (for 7^ = jl = 1 ) without 
interaction ||. Indeed, in the limit of large tK ^ 1 which 
corresponds to the unitary, pure ballistic case, our ex- 
pression for the current ( |l5| ) reduces to that derived in 
Rcf. H . This agreement is further supported by our nu- 
merical calculations carried out for tK = 100 (see Fig.l). 
Calculations of the current and conductance were also 
performed for tK = 5,3, and 2. The results are displayed 
in figure 1. For voltages in the sub-gap region, the con- 
ductance of an S(p-wave)KN junction is distinct from 
that of an S(s-wave)KN junction. Such difference shows 
that the ZES supports the formation of a Kondo sin- 
glet for the lower values tx = 5, 3 and 2 and effectively 
turn the junction to be more transparent, approaching 
the BTK limit §. 



Let us briefly summarize our results for the current 
and conductance of SKN junctions for s and p wave su- 
perconductors. In the limit of large tx the G—V curve is 
practically independent of tK and resembles that of a bal- 
listic junction, as indicated in the upper curve in Fig.l. 
For lower values of tK, junctions with s-wave supercon- 
ductor lead are driven away from the unitary limit and 
the underlying physics becomes much richer. It reflects 
the influence of both A and V on the Kondo resonance 
and on the actual value of the Kondo temperature. For 
tK = 5 the G—V curve noticeably deviates from that ob- 
tained in the non-interacting limit. For tK = 2 the com- 
petition between gap-related suppression of the Kondo 
effect and the effective transparency of the junction be- 
comes essential, leading to further decrease of the con- 
ductance. However, it is interesting to note that for and 
S(s-wave)KN junction at tK = 2 the conductance dis- 
plays a small peak at the gap edge. Its interpretation is 
that the Kondo correlations strongly compete with super- 
conductivity and influence the quasiparticle correlations 
when the energy exceeds the gap. Such an effect takes 
place even when tK < 1 (sec On the other hand, 

the S(p-wave)KN junction is less influenced to variations 
of tK- This may be explained by the fact that due to the 
presence of ZES, superconductivity plays a minor role in 
the formation of the Kondo resonance and BTK behavior 
persists for smaller values of tK- 
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SHOT-NOISE 



The shot-noise spectrum is usually dehned as the sym- 
metrized current-current correlation function M 



K{tx,t 2 ) = h[{I{tx)I{h)) - {If 



(17) 



where (...) denotes quantum averaging with the Hamil- 
tonian (Q) . The general expression for the zero frequency 
shot-noise power in junctions with one (or even two) su- 
perconducting leads (in MFSBA approximation) was ob- 
tained in our recent work fig] . Like in the case of SKS 
junctions, it is useful to write K — (K\ + K-^e 1 A/(8h) 
for which the expressions derived are, 



g*)(G R -G A )-g«G K }, (18) 




eV/A 



Fig. 2 The shot-noise power K (in units of e 2 A/h) as a 
function of V (in units of A/e) for an SKN junction. The 
parameters and notations are the same as in Fig. l.(tx 
grows from top to down. The curves with tx — WO coincide 
for both s and p wave superconductors. 



{Xt K f 



K 2 = -±lj^Tr{(G K ~gj 2 - 2r z ~gr z G A gG R 
[2gG R r z g K G K + (g R G R f - {G A g K r z f + h.c.p9) 



Expressions ( [l8| ) and ([II]) (supplemented by the self- 
consistency eqs. (|l^) and (M)) are then solved nu- 
merically for the same set of parameters T/T K = 200, 
t k = 100, 5, 3 and 2. The results for the shot-noise power 
spectrum K versus the applied voltage V are displayed in 
figure 2. These results are clearly correlated with those 
for the G — V curve and can be summarized as follows: In 
the limit tK 3> 1 the characteristics of shot-noise power 
spectrum for both s and p wave superconductors are con- 
sistent with those obtained for purely ballistic junctions 
which exhibit strong suppressing of the shot-noise power 
in the sub-gap region. At lower tK the physics is distinct. 
For tx = 5 the noise spectrum for s wave superconduct- 
ing lead still shows features typical for a junction with 
relatively high transparency, while the results for tK = 2 
are more similar to those for a low transparency junc- 
tion. Such dependence is explicitly exposed in the plot 
of the Fano factor versus the applied voltage (see Fig. 3). 
Though the Fano factor does not reach the maximum 
value of 2, it is strongly enhanced for the smaller value 
of tK = 2. For p — wave superconductor the shot-noise 
power (as function of voltage) reflects the same physics 
as in the conductance: ZES makes the Kondo resonance 
less vulnerable to the impact of superconductivity and 
the junction remains close to the unitary limit nearly in 
the whole range of values of tK considered here (see the 
solid curves on Fig. 2 and Fig. 3). 



1.0 - 




0.0 0.5 1.0 1.5 

eV/A 



Fig. 3 The Fano factor K/2eI as a function of V (in units of 
A/e) for an SKN junction. The parameters and notations 
are the same as in Fig. 1. 



JOSEPHSON CURRENT 

In this section we study an equilibrium property 
(Josephson effect) of an SKS junction in which both 
electrodes are p — wave superconductors (for comparison 
we also represent the results for s-wave superconducting 
leads). In the Kondo regime, the self-consistent equations 
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for the p — wave case(jl3|) and dll) now read 



2r ? 

— log — 

7T T K 



A 2 + /3(w)A 2 /wcos 2 | 



(l + a(w)) 2 w 2 



■/?H 2 A 2 cos 



2 i 



2r|w| 



2^ 



(1 + a(w)) 2 w 2 + e 2 + /3(w) 2 A 2 cos 



2 A 

2 



(20) 



(21) 



Here a(uS) = T^Jlu 2 + A 2 /iv 2 , /3(lj) = T/uj and S is the 
phase difference between two superconductors. Self con- 
sistent equations for the case of s — wave superconducting 
leads were derived in Ref. |2C| ] 

The Josephson current is obtained by variation of free 
energy with respect to the phase of the order parameter, 



Explicitly we then find, 



E 



2edF 
Ji~d6' 



(/3(w)A) 2 sin<5 



u 2 {\ + u(uj)) 2 + V + (/3(w)A) 2 cos : 



2 S 



(22) 



(23) 



The self-consistency equation and the expression for the 
Josephson current can easily be extended to the case 
of an anisotropic coupling '^l/r = T(l ±p). For this, 
one should replace cos 2 | — > (cos 2 | + p 2 sin 2 |) and 
/ oc sin S —y I oc (1 — p 2 ) sin<5. 

For a small anisotropy parameter p = 0.1 and for val- 
ues of tx =100, 5, 3, 2, we calculate the Josephson current 
as a function of the phase difference <5 at temperatures 
close to T = 0. The plot is depicted in Fig. 4. 



t«=100 




Fig. 4 Josephson current versus phase difference S at 
T — > 0. Dash and solid curves correspond to s and p-wave 
superconductors, correspondingly.The parameters are the 



same as in Fig.l. The value of tn decreases from top 
downward. 

In the s — wave case (dashed lines in figure 4) the cur- 
rent at tx =100 qualitatively corresponds to the unitary 
limit [M , while its amplitude is decreased at smaller val- 
ues of tx- Nevertheless, in this region of parameters, the 
Kondo effect strongly competes with superconductivity 
and therefore there are no traces of a 7r-junction. The 
mid-gap states for p—wave superconductors act as if they 
increase the effective normal region of the junction and 
therefore, the weight of the contribution to the Josephson 
current which is proportional to sind is reduced. When 
tx becomes smaller, the Josephson current for p — wave 
junction tends to be similar to that in an SIS junction. 




T/T c 

Fig. 5 Temperature dependence of the maximum Joseph- 
son current for s-wave (dashed curves) and p-wave(solid 
curves) cases when =100, 5, 3 and 2 (from the top 
to the bottom). The phase 8 of each curve is chosen 
so as to give the maximum value of Josephson current at 
T = in Fig. 4. The rest parameters are the same as in Fig. 1. 

By the same reasoning the temperature dependence 
(see Fig. 5) of the Josephson current is more similar to 
the usual SIS junction for the p — wave case, while for 
s — wave superconductors the current behavior is closer 
to the temperature dependence of the Josephson Current 
in SNS junctions. To be more precise, Fig. 5 displays 
the temperature dependence of the maximum Josephson 
current I(T) for s-wave (dashed lines) and p-wave (solid 
lines) cases. For p-wave junctions, I(T) assumes rela- 
tively large value near T = Tc while for s-wave junctions, 
it rapidly deceases as T increases. This difference stems 
from the question of whether the mid-gap state appears 
or not. The situation is similar to the one encountered in 
Josephson current through in an SIS system 27 . The 



difference becomes more prominent when the Kondo ef- 
fect is suppressed; I(T) drops more rapidly for s-wave 
junctions than for p-wave ones as tx decreases. 

In conclusion, we have analyzed an important physical 
problem involving strong correlations, the Kondo effect 
and superconductivity. These aspects can be combined 
in an SKN junction consisting of an Anderson impurity 
(in the Kondo regime). We have developed a theoreti- 
cal framework by which it is possible to investigate an 
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interplay between Andreev reflections and the formation 
of Kondo resonance in the Kondo regime T < A < Tjc . 
In this limit we calculated Josephson Current. We have 
also investigated non-equilibrium aspects and elucidated 
the nonlinear G—V characteristics. Moreover, we calcu- 
lated the shot-noise power spectrum of SKN junctions 
at voltages eV < A. It is found that at sufficiently large 
tx the Kondo resonance effectively turns the junction 
behavior to be similar to that of highly transparent non- 
interacting weak links for both s and p — wave super- 
conductors. However, when the ratio of the Kondo tem- 
perature to the superconducting gap becomes smaller the 
behavior of these two types of junctions is quite different: 
the Kondo resonance persists much more effectively for 
junctions with p — wave leads than it does for s — wave 
leads. 
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